3D view of transient horizontal magnetic fields in the photosphere 
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ABSTRACT 



We infer the 3D magnetic structure of a transient horizontal magnetic field 
(THMF) during its evolution through the photosphere using SIRGAUS inversion 
code. The SIRGAUS code is a modified version of SIR (Stokes Inversion based 
on Response function), and allows for retrieval of information on the magnetic 
and thermodynamic parameters of the flux tube embedded in the atmosphere 
from the observed Stokes profiles. Spectro-polarimetric observations of the quiet 
Sun at the disk center were performed with the Solar Optical Telescope (SOT) 
on board Hinode with Fe I 630.2 nm lines. Using repetitive scans with a cadence 
of 130 s, we first detect the horizontal field that appears inside a granule, near 
its edge. On the second scan, vertical fields with positive and negative polarities 
appear at both ends of the horizontal field. Then, the horizontal field disappears 
leaving the bipolar vertical magnetic fields. The results from the inversion of the 
Stokes spectra clearly point to the existence of a flux tube with magnetic field 
strength of ~ 400 G rising through the line forming layer of the Fe I 630.2 nm 
lines. The flux tube is located at around logrsoo ~ at At=0 s and around 
logr 50 o ~ —1.7 at At=130 s. At At=260 s the horizontal part is already above 
the line forming region of the analyzed lines. The observed Doppler velocity 
is maximally 3 km s _1 , consistent with the upward motion of the structure as 
retrieved from the SIRGAUS code. The vertical size of the tube is smaller than 
the thickness of the line forming layer. The THMF has a clear fi-shaped-loop 
structure with the apex located near the edge of a granular cell. The magnetic 
flux carried by this THMF is estimated to be 3.1 x 10 17 Mx. 

Subject headings: Sun: magnetic fields — Sun: photosphere 
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Introduction 



The SOT spectropolarimeter 



Ichimoto et al 



2008; 



Shimizu et al. 



f SP. iTsuneta et al.l 12008a 



20081 ) on board Hinode (IKosugi et al. 



Suematsu et al.l 12008c 



20071 ) revealed 



that there are ubiquitous horizontally i nclined magnetic fields with strengths o : 



hundred G in t he in terne twork regions ( Orozco Suarez et al. 



Centeno et al. 



(120071 ) and 



Ishikawa et al. 



2007a ; 



Lites et al. 



a few 



2008|). 



(120081 ) reported the temporal evolution of such 
horizontal magnetic fields both in the quiet Sun and in plage regions. The evolution 
commonly starts with the appearance of linear polarization signals (horizontal fields) inside 
granular cells, followed by the appearance of the opposite-polarity circular polarization 
signals (vertical fields) at both ends of the horizontal component. The time evolution 
indicates that the horizontally inclined magnetic fields are part of a small magnetic loop. 
The sizes of the horizontal magnetic fields are smaller than the size of the granule where 
they appear. Statistical analyses of a large number of the horizontal fields indicate that 
the life time of the horizont al magnetic fields ranges fro m 1 to 10 minutes and the mean 



lifetime is about 4 minutes (jlshikawa fc Tsuneta 



2009a 



bj). Their occurrence rate is high, 



and more than appro ximately 10% of the granules have embedd ed horizontal magnetic fields 



( jlshikawa et al. 



20081 ) . The movie of linear polarization signals (jlshikawa fc Tsunetal l2009al ) 



shows that these horizontal magnetic fields are highly frequent and transient, and the term 



'transient horizontal magnet ic fields" (THM 



emer gence in the quiet Sun ( jLites et al 



1996 



s) is used for them. Earlier reports of flux 



DePontieu 



2002 



Martinez Gonzalez et al. 



20071 ) appear to be examples of THMFs. 



The properties of the THMFs described above are the same in the quiet Sun and 
the plage region. There are granular-sized horizon tal fields in the polar regions, and 



these magnetic fields may be the same as THMFs ( ITsuneta et al. 



2008b 



Ito et al. 



20091 ). 



However, there seem to be a few distinct patterns in their disappearance. For example, 
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lane, while 



Ishikawa et al.l (120081) rep o rt an example of a THMF disappearing in the intergranular 



Centeno et al. 



This is confirmed by 



(200 71) show a T 



TMF which disappears inside a bright granule. 



Ishikawa fe Tsunetal (l2009bl ) who report that normalized continuum 



intensities where THMFs appear are above ~1.0, corresponding to the bright granular 
regions. However, the normalized continuum intensities where THMFs disappear are 
both higher and lower than ~1.0, i.e. THMFs do not necessarily reach the intergranular 
lanes. There also appears to be two types of THMFs; with or without apparent bipolar 
vertical fields after their disappearance. This variety of THMF behavior suggests a few 
different scenarios with regard to their disappearance: (1) THMFs somehow submerge 
with the downward convective motion; (2) They become fragmentary due to the convective 
flow, and their Stokes signals are below the SP detection level; (3) THMFs go through 
the line formation layer of the Fe I 630.2 nm lines and reach higher atmospheric layers 



( IMartmez Gonzalez fe Bellot Rubio 



20091 ). 



Given the ubiquity of THMFs in the solar photosphere, the question of how 
they appear and disappear in such a short time becomes one of the important issues 
in sol ar magnetohydrodynami c s (MHD). A series of sophisticated MHD simulations 



Vogler fe Schiissler 



Steiner et al 



2007 



Abbott 



2007 



Schiissler fe Vogler 



2008 



Isobe et al. 



2008 



20081 ) have been carried out with varieties of initial and boundary conditions. 



A detailed comparison of the observed properties of THMFs with the numerical simulation 
results has just begun. It is thus desirable to clarify the structure and evolution of the 
enigmatic small-scale magnetic fields in the photosphere, and possibly in the chromosphere. 



Spectropolarimetric data obtained with the SOT have been anal yzed thus 



with inversion codes based on a Milne-Eddington atm osphere (e.g. 



Orozco Suarez et al 



2007a 



' ar ma inly 



Lites et al. 



2008 



Ishikawa fe Tsunetal l2009al ). Milne-Eddington inversion 



provides us with various physical parameters (magnetic field vector and Doppler velocity, 
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etc.), which have to be constant throughout the atmosphere (except for the source function), 
and synthesize only symmetric Stokes profiles. However, the observed Stokes profiles of the 
horizontal magnetic fields usually have area-asymmetries. Such asymmetric profiles suggest 
the presence of velocity and magnetic field gradients a long the line of sight (LOS) (e.g., 



Solanki fc Pahlke 



1988 



Sanchez Almeida Lites 



19921 ). Thus using the Milne-Eddington 



inversion, physically important information about the small-scale horizontal magnetic fields 
may be lost. 

We are usin g an inversion code based on the SIR code (Stokes Inversion based on Re - 



spon s e function, 



1996 



Ruiz Cobo fc del Toro Iniesta 



Del Toro Iniesta 



1992 



1994 



Del Toro Iniesta &: Ruiz Cobo 



20031 ). This code allows for the change of atmospheric pa- 



rameters along the LOS and thus fits the observed asymmetric Stokes profiles. 



Martinez Gonzalez et al 



( 120071 ) used the SIR inversion code to analyze a snapshot of the 
magnetic loops using the spectropolarimetric data of 1.56 fim taken with the Tenerife 
Infrared Polarimeter. However, the used settings of the inversion assumed the plasma 
parameters to be constant with heigh t in the atmosphere, except in magnetic field 



inclination and temperature. Recently, 



Gomory et al. 



(120091 ) studied the temporal evolution 



of a similar event applying SIR inversions. However, they kept all magnetic parameters 
constant with height. Thus, these analyses resemble a Milne-Eddington inversion with 
a realistic temperature distribution. We fully exploit the inversion code without such 
constraints to analyze the observed Stokes profiles asymmetries and with a better spatial 
resolution. This allows us to follow the temporal evolution of the stratified atmosphere. 



by 



In this paper , we st udy a THMF event whose properties are similar to those reported 



Centeno et al 



(20ol 



code (jBellot Rubio 



. We carry out the full Stokes inversion with the SIRGAUS 



20031 ). which is a modified version of SIR. The SIRGAUS code treats 



explicitly a magnetic flux tube embedded in the atmosphere with a Gaussian function, and 
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is expected to be more suitable to obtain vertical (LOS) structure of the THMF. Indeed, we 
successfully identify an isolated flux tube in the photosphere, and track the emergence of 
the f2-shaped flux tube through the photosphere. We describe statistical significance of our 
result in detail. With the clear identification of the upward moving flux tube, we are able 
to obtain various physical parameters that characterize the flux tube. Such information on 
the flux tube becomes available for the first time with the use of the SIRGAUS code. We 
finally discuss the magnetohydrodynamic properties of the THMF. 



2. Observations and wavelength calibration 

A quiet solar region of 4" x 76" located near the Sun center was observed with a 
cadence of 130 s using the SOT/SP on March 10, 2007. We obtained the Stokes spectra 
of the Fe I 630.2 nm lines, w hose line formation region is between logrsoo ~ and —2.5 



(jCabrera Solana et al. 



20051 ). Each slit position has an exposure time of 4.8 s, and the 
scanning step is / .'15. The noise level in the Stokes Q and U signal is 1.2 x 10 _3 / c and 
that in the Stokes V is 1.1 x 10~ 3 / c , where I c is the continuum intensity. The dark current, 
the flat field corrections, and the polarimetric calibration were performed using "sp_prep" 
software. 

We extensively analyze the event shown in Figure [TJ The horizontal magnetic field 
appears at the edge of a granule without any vertical magnetic fields (At = s). At 
At = 130 s (second frame), the size of the horizontal magnetic field becomes larger and 
the bipolar vertical fields appear at both ends of the horizontal field. At the third frame 
(At = 260 s), the horizontal field disappears, leaving an isolated pair of bipolar vertical 
fields. The separation of these two patches increases with time. The lifetime of the 
horizontal magnetic field is about 4 min, and this event has properties in common with 
those of the THMFs as reported by 



Ishikawa et al. 



fl2008h 
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We analyze only the pixels with the Stokes Q or U amplitude 10 times larger than their 
noise level. Such pixels are shown with the red contours in Figure [TJ Pixels shown with the 
green contours are dominated by the Stokes V signal, and correspond to the footpoints of 
the THMF. These pixels are also analyzed. 

The Doppler velocity information is crucial for understanding the dynamics of the flux 
tube. For this purpose, accurate in-orbit calibration of the SOT/SP wavelength scale and 
offset is done to obt ain the absolute velocity scale. We here calibrate our data using the 

i — — inn 

solar flux atlas (IBeckers et al.lll976l ). The atlas contains a high resolution spectra of the 
solar irradiance from 380 nm to 700 nm and has the wavelength information. 

First, we obtain the mean Stokes / profile of the quiet Sun by averaging the Stokes 
I profiles of pixels without the polarization signal. Second, we compare the mean Stokes 
I profile with the solar flux atlas and derive the wavelength sampling of 2.14 pm and 
the wavelength value at the start pixel of 630.089 nm. These two parameters give us 
the maximum correlation coefficient between the mean Stokes I profile obtained with 
the SOT/SP and the solar flux atlas (see Fig. [2]). However, the wavelength shift of 



the flux atlas is not corrected 



or the gravitational redshift and convective blueshift 



(lAUende Prieto fe Garcia Lopez 



19981 ). It turns out that the two Fe I lines have the 
redshifts of 0.88 pm and 0.63 pm (420 m s _1 and 300 m s _1 ) with respect to the laboratory 
centers, respectively. These different Doppler shifts come from the difference in the 
convective blueshifts between the two lines. Thus, we define the mean value of the two 
redshifts for the two Fe I lines to be 360 m s _1 as the zero point of the velocity scale. With 
this process, we have corrected the effect of the gravitational redshift and the convective 
blueshift. The Doppler velocities in this paper are now defined such that if there is no 
motion on the Sun, this will give the zero velocity. 
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SIRGAUS Inversion 



We perform th e inversion of the observed Stokes profiles using the SIRGAUS code 



f lBellot Rubio 



20031 ). This inversion code synthesizes the Stokes profiles emerging from a 
model atmosphere and then compares them to the observed ones. Using the least-square 
Levenberg-Marquardt algorithm, the initial model is modified until the difference between 
the observed and the synthetic profiles is minimized. The model consists of a simple 
atmosphere that represents the background and a Gaussian-shaped bump superposed on 
some of the atmospheric parameters representing an embedded flux tube. Note that we are 
considering a one-component model atmosphere. The THMF analyzed here is spread out 
over several pixels. Therefore, we assume that the magnetic filling factor, /, is one. 

Recent Milne-Eddington inversions o f THMFs have provided an average filling factor 



of about 0.2 ( jlshikawa fc Tsuneta 



2009aJ). The observed magnetic structures we want to 
analyze cover more than a spatial resolution element of the S P data. Thus, if we neglect 
stray-light contamination effects like those due to diffraction ( Orozco Suarez et al. 2007bl ) 
or scattered light, we would expect / to be around unity. The obtained small filling factors 
might be just a net effect of the used technique that cannot account for a flux tube fully 
covering the resolution element, but only partially occupying the line forming region of the 
Fe I lines at 630.2 nm. 

Thus far, this code has not been widely used. Examp l es of its application to 



Jurcak et al. 



(120071 ) where they analyze 



spectropolarimetric data can be found in, e.g 
Stokes profiles observed in Sunspot penumbrae. As explained before, in SIRGAUS the 
distribution of the atmospheric parameters, x, consists of a background atmosphere and of 
a Gaussian bump, and can be written as 

(V ~ Gpos) 2 



x(n) = x (rj) + Axexp 



C 2 

width 



1 
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where x b stands for the background plasma parameters; Ax represents the amplitude of the 
Gaussian bump; 77 is the logarithm of the optical depth, r 500 ; and G wi dth an d G pos are the 
width and the position of the Gaussian bump in rj units. The value of the amplitude, Ax, 
is different for each of the physical quantities. The Gaussian bump position and width are 
the same for all the physical quantities. 

The SIRGAUS code tries to find the best model atmosphere that reproduces the 
observed profiles in a space of 14 free parameters: The temperature, T, of the background 
atmosphere is allowed to change at three nodeJl]. The other free parameters, x b , 
characterizing the background atmosphere, say the microturbulence, v b mic , the LOS velocity, 
v losi an d the field strength, B b , its inclination, j b , and azimuth, (fi b , are assumed to be 
constant with height. The rest of the free parameters are the amplitudes of the bumps, 
Avlos, AB, A7, and A0, and their width, G W idth, and position, G pos . Hereafter and 
for simplicity, since the Gaussian bump describes the plasma parameters of a flux tube 
embedded in the line forming region of the observed spectral lines, the Gaussian bump will 
be referred to as flux tube. 

Throughout this paper, we are displaying the resulting distributions of plasma 
parameters as a function of the logarithm of the continuum optical depth at 500 nm; 
logT5oo, hereafter referred to as logr. However, the optical depths will be converted to 
kilometers to specific purposes (see e.g., Sect. I6.4p . To this end we solve the following 
equation: r = L updz. With the boundary condition r = 1 at z = km, k is the continuum 
absorption coefficient per unit mass and p the gas density. The absorption coefficient 
depends on different plasma parameters and on the chemical composition. The results 
from the SIRGAUS inversion are used to evaluate the previous equation at each pixel. For 
its integration we use the program geometrical. x, provided with the SIR code. The same 



x See e.g., 



Ruiz Cobo fc del Toro Iniestal (119921 ) 
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program also provides us with the values of the gas density that will be used later in this 
paper (see Sect. I7.2p . 



In order to clearly show the evolution of the horizontal magnetic field between the 
two scans, we select two pixels, one in each scan, and present their detailed analysis. The 
selected pixels are marked by yellow crosses in Figure [TJ The corresponding observed Stokes 
profiles are shown in Figure [31 To invert the Stokes profiles of all pixels we used several 
initial atmospheric models. For the pixel taken from the first scan (hereafter pixel^l), 
all the initial models converged to a solution where the flux tube is located low in the 
atmosphere. On the other hand, the solutions for the pixel from the second scan (pixel#2) 
show that the position of the flux tube is located higher in the photosphere. 

To show the reliability of the obtained positions and widths of the flux tube, we 
proceeded as follows: First we fixed all the plasma parameters obtained by the SIRGAUS 
inversion. Then we evaluated a merit function for a wide range of G W idth and G pos 
values. For this purpose, we introduced a x 2 function that does not take into account the 
contribution of the Stokes / profile: 



where n is the number of the wavelength points; Ik with k — 1,2, 3 indicates the Stokes 
Q, U, and ^profiles; I% bs and l s k yn stands for the observed and the synthetic profiles, 
respectively; and a Q is the standard deviation of the Stokes Q, U, and V profiles estimated 
using pixels without any apparent polarization signal. It represents the Poisson noise of 
incident photons and the value is ~ 1.1 x 10~ 3 / c (See Sect. [2]). 

Hereafter, we use the y 2 function as defined in Eqn. (T5]). The SIRGAUS code minimizes 



4. Emergence of horizontal magnetic fields 




3 



(2) 



fc=i i=i 
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a x 2 that takes into account all four Stokes profiles. The synthetic Stokes profiles shown in 
Figure [3] are example fits to the observarions. As can be seen in the lower panels (residuals), 
the fits are very good for the Q, U, and V profiles and reasonable for the / profile. The 
largest discrepancies are obtained in the wings of the Stokes I profile with a maximum of 



taking into account the hydrogen collisional broadening (jBellot Rubio 



s of Stokes I 


may fc 


(IBellot Rubio 


2009) 



should be noted that SIRGAUS retrieves the information about the atmospheric structure 
along the LOS mainly from the asymmetries of the Stokes Q, U, and V profiles. For this 
reason we use a \ 2 based only on the Stokes Q, U, and V profiles. This x 2 definition will 
also help in clarifying the properties of the flux tube. 

Figure H] shows the dependence of the \ 2 values on the position and width of the flux 
tube for pixel^l and pixel#2 (left and right panels, respectively). The minimum values 
(Xmin) are marked by the crosses and correspond to the synthetic profiles shown in Figure [3] 
and to the models shown by the gray lines in Figure |5j To construct Figure HI we fix 12 of 
the 14 free parameters and evaluate the \ 2 value using Eqn. (T5]) and varying the G pos and 
Gwidth- The solid curves in Figure 0] are defined by constant values of Ax 2 = X 2 ~ Xmin- 
The region outlined by the Ax 2 curves defines the level of confidence: for instance, there is 
a 68.3% (la) chance for the true parameters of the position and width of the flux tube to 
fall within the region defined by the yellow contour (A% 2 = 13.7). Likewise, there is 95.4% 
(2a) and 99.7% (3a) chance that the true parameters fall inside the red contours with 
Ax 2 = 21.30 and 29.8, respectively. As explained before to make Figure HJ we have varied 
2 parameters (the positions and widths of the flux tube) and remaining 12 parameters 
are fixed. In this case, Ax 2 is dis tributed as a chi- square distribution with 14 — 2 = 12 



degrees of freedom (TheoremD in 



Press et al. 



20071 ). The A X 2 =13.7, 21.3, and 29.8 used 



above come from the chi-square distribution function for 12 degrees of freedom and for the 
confidence level of 68.27%, 95.45%, and 99.73%, respectively. 
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In pixel^l (At = s), the \ 2 is minimum when the position of the flux tube is about 
logr ~ —0.15 and its width is about 1.1. However, these parameters are not determined 
uniquely since the la confidence level is elongated. The positions of the flux tube extend 
from logr ~ 0.1 to —0.7 with G W idth changing accordingly from 1.4 to 0.5. A thick flux tube 
located deeper in the photosphere and a thin flux tube located higher in the atmosphere 
provide equally good fits to the observed spectra. From this results we can ascertain that 
the center of the flux tube is not located higher than logr ~ —0.7. This optical depth 
corresponds to a geometrical height of 91 km. 

The lcr confidence level is less elongated for pixel#2 (At = 130 s), but again, it is 
difficult to specify the exact width and position of the flux tube. The minimum x 2 is 
obtained when the position of the flux tube is about logr ~ —1.8 and its width about 
1.1. The position of the flux tube is definitely higher than logr ~ —1.5, i.e., 213 km. 
In conclusion, the \ 2 maps clearly show with sufficient statistical significance that the 
magnetic flux tube is located higher in the atmosphere at At = 130 s than at At = s. 
Therefore, the flux tube moves upwards with time. 

The distributions of plasma parameters as a function of the optical depth are shown 
in Figure [5j The gray line stands for the atmospheric model at Xmin- The black areas 
represent the variation within the la confidence region. We first notice that the flux tube 
does not occupy the whole line forming region. For pixel^l, the center of the flux tube is 
around logr ~ —0.5 and we are observing the upper boundary of the flux tube. On the 
other hands, for pixel#2, the flux tube is located at around logr ~ —1.7 and the lower 
boundary of the flux tube is observed. For both pixels, the background plasma is not 
magnetized even though we allow the field strength to change in the inversion. Thus, the 
flux tube is located in an essentially field-free atmosphere. 

The plasma parameters such as the magnetic field strength significantly change as a 
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function of the optical depth in the line forming region due to the presence of the flux 
tube (Fig. [S]). The large gradients are probably responsible for the observed asymmetries 
of the Stokes profiles. We note that the variation in the plasma parameters (as represented 
by the black areas in Fig. [5]) at the same optical depth is relatively small, regardless of 
the considerable differences in the widths and the positions of the flux tube, essentially 
within the line forming regions (the region delimited by the two dashed lines in Fig. For 
pixel^l, the values above logr ~ are similar: the broader the flux tube becomes, the 
deeper is located to maintain the similarity in the distribution. Likewise, pixel#2 has a 
similar situation. The values of the plasma parameters are similar at the same optical depth 
below logr ~ —2.0, regardless of the widths and the positions of the flux tube. Broader 
flux tubes are positioned higher in the atmosphere for the same reason. 

The rising motion of the flux tube is confirmed by the distributions of LOS velocity. 
The LOS velocity distributions shown in Figure |5] indicate that the flux tube is slightly 
blueshifted below logr = —0.7 for pixel^l. The velocity of the flux tube shows an upward 
motion with respect to the background atmosphere at logr ~ —1.5. The relative velocity 
of the flux tube is ~ 1.3 km s -1 . For pixel#2, the flux tube is blueshifted, showing LOS 
velocities up to 1.6 km s _1 . The relative velocity with respect to the background is larger 
than for pixel^l. 

As can be seen in Figure the inclination angle of the flux tube is ~ 90°, so that the 
field lines are perpendicular to the LOS for both pixels. The azimuth angle, measured from 
the east- west direction, is about 160°; indicating that the field lines are nearly parallel to 
the line connecting the negative and positive Stokes V patches of the second scan. The 
magnetic field strengths are around 400 G in both cases. 
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5. Three dimensional structure and time evolution 

We have analyzed those pixels outlined by red and green contours (See Fig. [I]), following 
the approach we applied to pixel^l and pixel#2 in the previous section. Figures [6] and [7] 
show the corresponding \ 2 maps for these pixels. Pixel^l and pixel#2 correspond to pixel 
1-6 in Figure |6]and pixel 2-14 in Figure [7J respectively. 

In the first scan, the results are comparable to those found for pixel^l. The la 
confidence regions are elongated and located low in the atmosphere. However, this is not 
the case for pixels 1-1 and 1-4. In these pixels, although the fits are as good as those in the 
neighboring pixels, the location of the flux tube is different. 

In the center of the THMF (pixels 1-6 and 1-2), the la confidence region is located 
between logr = and logr = —0.7. There the position of the flux tube is the highest. In 
the surrounding pixels (1-3, 1-5, 1-7, 1-9, 1-10, 1-11) the la confidence region is located 
below logr ~ 0. At the one of edges of the THMF (pixels 1-8 and 1-12), the location of 
the flux tube is even lower. In summary, if we go from the center of the THMF to any 
of its edges, the flux tube sinks into the photosphere. Notice that for most pixels the 
lower boundary of the flux tube is below the continuum formation level. This results is 
compatible with the presence of a low laying flux tube with an Q shape. Only the apex 
of the flux tube is detected in linear polarization (i.e., horizontal fields) and therefore, the 
footpoints of the loop are not seen. 

In the second scan, in the center of the THMF, the flux tube is located at around 
logr = —1.6. When we go from the center of the THMF towards any of its footpoints (e.g., 
from pixels 2-15 to 2-19), the flux tube is located deeper in the atmosphere. Again, this is 
consistent with a flux tube with a Q shape. The retrieved shape and position of the flux 
tube from the first and second scan confirms our observation of a rising Q loop. 
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6. Properties of transient horizontal magnetic field 

6.1. Magnetic field topology 

The lower left part of Figures [6] and [7J shows the azimuth angles of the magnetic field 
lines retrieved with the SIRGAUS code at logr = —1.0. For both scans, the magnetic field 
is aligned with the axis of the flux tube, i.e., with the line connecting the opposite polarity 
patches. The averaged azimuth angles around the center of the THMF are 164° and 162° 
for the first and second scan, respectively. Since the azimuth is displayed at the same 
optical depth for both scans, it means that the direction of magnetic field is the same in 
the upper part of the flux tube observed at At = s as in the lower part of the flux tube 
observed at At =130 s. This indicates that the analyzed THMF does not show any helical 
structure, so that essentially, the flux tube has a toroidal component only. 

The pixels in which the azimuthal directions are shown with green lines in Figures |6] 
and [7J are selected for plotting the magnetic field strength cross-section of the flux tube, 
along its axis (top panels in Fig. Left and right panels correspond to the first and second 
scans. For both scans, the field strength distributions clearly show that the flux tube is 
located deeper in the atmosphere, going from the center to the edge of the THMF. In the 
first scan, we detect only the apex of the flux tube, while in the second scan the flux tube 
is already located high in the line forming region. These results clearly indicate a rising 
fi-shaped flux tube. The magnetic field strength of the flux tube is smaller than 400 G. We 
do not see any change in magnetic field strength between the first and second scan. 

The white lines plotted over the distributions of the magnetic field strength in Figure [8] 
mark its inclination. In the first scan only the horizontal component of the magnetic field 
can be detected, 7 ~ 90°, (there is no Stokes V signal observed during this scan). The 
inclination at the central part of the THMF is also around 90° in the second scan. However, 
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near the region where Stokes V signal is detected, the retrieved inclinations indicate that 
the field lines are more vertical. The retrieved inclinations nicely delineate the shape of the 
flux tube. 



6.2. LOS velocity topology 

The distributions of LOS velocity along the axis of the flux tube can be seen in Figure El 
(middle panels). Since the plasma velocity depends on whether the pixel is located in the 
granular /intergranular region, the normalized continuum intensities are also displayed in 
the bottom panels of Figure 

Pixels 1-4 and 1-5 have continuum intensities lower than 1.0, and correspond to the 
edge of the granule (see also Fig. [1]). In these pixels, the flux tube show downward motion. 
In the same pixel, the background atmosphere shows larger downflows. From pixel 1-10 
to pixel 1-12, the continuum intensity is above 1.0 and becomes larger. These pixels 
correspond to the transition region from the edge of the granular cell to the interior of the 
granule. There, the flux tube exhibits a strong upward motion, with velocities reaching 
2.7 km s -1 . In contrast, the background plasma above the flux tube is slightly redshifted 
with LOS velocities smaller than 1 km s _1 . Note that t he absolute Doppler velocities 
with respect to the solar surface are shown in Figure [8j iBorrero &: Bellot Rubiol ( 120021 ) 
obtained a two-component model of the solar photosphere from the inversion of spatially 
and temporally averaged solar intensity profiles, and they found upflows of ~1.1 km s -1 
in granular cells at around logr ~ —0.3 (Note that they only correct the gravitational 
redshift). Thus, the upward velocities of the flux tube we find in our analysis are larger 
than that for typical granular cells. 

In the second scan, the flux tube is located higher in the atmosphere. At the apex of 
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the loop (pixel 2-13), the continuum intensity is slightly below 1.0. This pixel appears to be 
located at the boundary between the granule and the intergranular lane. The pixels around 
the apex of the flux tube show upward velocities reaching ~ 3 km s -1 , clearly indicating 
the rising motion of the flux tube. 

Next, we examine the more vertical portion of the flux tube located in the low 
atmosphere. The region of negative polarity (pixels 2-17, 2-24, 2-25, and 2-26) shows the 
continuum intensities larger than 1.0, corresponding to the inner part of the granule. There, 
the flux tube has strong upward velocities ranging from ~ 2 km s _1 to ~ 4 km s _1 . These 
velocity values are larger than that of the typical upflows observed in granules. 

In the portion of the flux tube with positive polarity (pixels 2-3 and 2-4), the 
continuum intensity is lower than 1.0. These pixels are located in an intergranular lane. 
The flux tube at these positions has a downward motion ranging from ~ 2.3 km s _1 to 
~ 4.6 km s _1 at —0.5 < logr < 0. Typical downflow velocities in the intergranular lanes 
i n the lower photosphere (—0 .5 < logr < 0) are between 2.2 km s -1 and 3.7 km s -1 



(IBorrero fc Bellot Rubio 



2002) comparable to what we measured. Therefore, the downward 



motion of the flux tube may be driven by the convective motion so that the flux tube may 
be forced to be submerged with the downward flow field. 

As we have seen, the flux tube is located at around logr = in the first scan (At = s), 
while in the second scan (At = 130 s), the center of the flux tube is located at around 
logr = —1.7. This clearly indicates that the flux tube goes upward and reaches higher 
photospheric layers with time. The geometrical distance between logr = and logr = —1.7 
is around 244 km. Given the time difference of 130 s between these two scans, the rise 
velocity of the flux tube is 1.9 km s _1 . This is consistent with the inferred upward motion 
of 2 km s _1 for the flux tube (See middle panels of Fig. 
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6.3. Cross sectional shape of flux tube 

To estimate the width of the flux tube, we need to find a pixel, where the la confidence 
regions in the x 2 maps (marked by the yellow contour in Figures [6] and [7]) are not elongated 
and where we observe both the upper and lower boundary of the flux tube. In the second 
scan, the center of the analyzed THMF is already too high for our purpose. We select pixel 
2-16. This pixel is located at the center of the THMF, and its x 2 icT confidence levels are 
not spanning over a wide range of G wi dth and G pos . The upper and the lower boundary are 
within the line forming region, as well (See the top right panel from Fig. |Sj). 

For pixel 2-16 the position of the flux tube is located at logr = —1.4 and the width 
corresponds to 0.8 in logr units as given by the Xmin- We define the radius of the flux tube 
along the LOS as the FWHM of the Gaussian function that represents it. In this case the 
vertical extension of the flux tube is around 190 km. We estimate the horizontal size to 
be approximately 360 km, since the THMF occupies three pixels (2-9, 2-16 and 2-22) at 
around the location of pixel 2-16. It should be kept in mind that the real horizontal size 
might be larger since a weaker Stokes Q and U signal can be observed surrounding the 
THMF. The discrepancy between the vertical and horizontal size of the flux tube indicates 
that the flux tube is squeezed in the vertical direction. 

6.4. Total magnetic flux carried by THMF 

In this section we estimate the magnetic flux carried by the THMF, since we know 
the vertical (along the LOS) and lateral cross section of the flux tubes as shown in the 
previous sections. Let B(z) be the magnetic field vector with inclination 7(2), z being the 
geometrical height, and the SOT/SP pixel size I' of about 120 km. The per-pixel magnetic 
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flux ifh through the surface S (see Fig. [9]) is given by 



B-dS 



2i(logr=-2.5) 



B{z) cos(90° -7(z))(l'dz). 



(3) 



'2 (logr=0) 

The integration range is set to cover the line forming region of the Fe I lines. To estimate 
the total magnetic flux carried by the THMF we sum up the (ph values for pixels of 
2-1, 2-7, 2-14, and 2-20 from the second scan. The total magnetic flux &h turns o ut to 



Jinetal 



(J2009|), who 



be 3.1 x 10 17 Mx. This value is in agreement with that derived by 
assumed that the vertical extent of the flux tubes is the same as that of the scale height of 
100 km. In our analysis we resolve the vertical structure of the magnetic flux tube using 
SIRGAUS. This has allowed us to derive the magnetic flux more directly. Note that the 
total magnetic flux may be different depending on the polarization threshold used in the 
analysis, so that our magnetic flux estimate represents a lower limit. 



7. Summary and discussions 

7.1. 3D magnetic structure and evolution 

We have analyzed a magnetic phenomena in which horizontal magnetic field patch 
appears first at the edge of a granular cell at At = s. Later at At — 130 s, the size of the 
patch harbouring horizontal fields becomes larger in size, and vertical magnetic fields with 
opposite polarity appear at both ends of the magnetic patch. Then (at At = 260 s), the 
horizontal magnetic field disappears, while the vertical magnetic fields remain and separate. 
The whole process takes about 4 min. This event corresponds to a typical example of the 
so called transient horizontal magnetic field (THMF). 

We have successfully applied the SIRGAUS inversion code to this event and analyzed 
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its 3D magnetic structure. The results show that the flux tube occupies only a fraction of 
the line formation region of the Fe I lines at 630.2 nm and that it has an upward movement. 
The magnetic flux tube rises from logr ~ (0 km) to logr ~ —1.7 (244 km) between the 
first and the second scan that are separated in time by 130 s. The Doppler velocity of the 
flux tube consistently indicates the upward motion. The magnetic field strength of the flux 
tube is about 400 G. The azimuth direction of the magnetic field is aligned with the axis of 
the flux tube, and the flux tube does not show any helical structure. 

So far, to determine whether the observed magnetic loops correspond to the rising 
fl loops or to submerging U loops, Dopp ler velocity measurement have been used (e.g., 



Martinez Gonzalez fe Bellot Rubio 



20091 ). Our analysis with SIRGAUS presented here 
shows that the flux tube with shallow but clear Q shape rises. We observe only a small part 
of the shallow Q loop above the continuum formation layer and this might be just a part of 
a larger-scale magnetic structure located in the bulk of the convection zone. 



7.2. Force balance of flux tube 

The field strength of the flux tube is ~400 G in the first scan (At = s). Taking the 
velocity of the convective flow to be v = 2 x 10 5 cm s _1 and the mean density at the base of 
the photosphere (logr ~ 0) to be p =2.7x 10~ 7 g cm -3 , we find that the equipartition field 
strength is B eq = (4m , p) 1 / 2 v = 370 G. Therefore, the inferred field strengths are comparable 
to B eq . This is consistent with the fact that the flux tube seems to maintain its integrity 
during its rise from the lower photosphere upwards. 

In this section, we compare the buoyancy force, Ff, = (p e — pi)g, with the magnetic 
tension force, Ft = B 2 /AtxL, acting on the flux tube at At =130 s. Here pi and p e are the 
plasma density inside and outside the flux tube respectively, g is the gravity acceleration of 
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2.7 x 10 4 cm s -2 , and L is the curvature radius of the Q loop. Here we use all the plasma 
parameters (temperature, field strength, and the density) at the height where the flux tube 
is located and the 3D geometrical scale of the flux tube (the length and the cross section). 
From the separation of the footpoints in the second scan (which is about 1320 km) and 
from the averaged magnetic field inclination in the footpoints (about 47°), we estimate L to 
be 1250 km. The comparison is made at the apex of the flux tube which is located above 
logr ~ —1.0. There the temperature T is ~4950 K (corresponding to the averaged value 
from logr = —1.0 to —2.0 at pixel 2-14). Assuming that the temperature is the same inside 
and outside of the flux tube, F b can be written as F b = |^ • jj^Sf, where m p is the proton 
mass, and ks the Boltzmann constant. The ratio between Ft, and Ft is then = ~ 2. 
This clearly indicates that the buoyancy force, Fb, participates in the emergence process of 
the magnetic flux tube. 

If we suppose that the stationary state is reached, then the buoyancy force, the 
magnetic tension force, and the resultant drag force associated with the rising motion are 
balanced: 

B 2 B 2 Tn q B 2 

(p e ~ Pi) 9 ■ al - -r-r • al - 0.5C d p e u 2 l = — ■ ■ al - — — ■ al - 0.5C d p e u 2 l 

4ttL 57T 2kb± 47tL 

= 0. (4) 



1979|), and u the 



With Cd being the aerodynamic coefficient of the order of unity (jParkerl 
relative velocity of the flux tube with respect to the velocity of the surrounding plasma, a 
and I are the vertical extent and the horizontal size of the flux tube, respectively. We take 
p e to be 7.78 x 10~ 8 g cm -3 (corresponding to the mean density from logr = —1.0 to —2.0 
at pixel 2-14). The relative velocity u obtained from Eqn. (j3J) is then 2.3 km s" 1 , which is 
consistent with the observed relative Doppler velocity associated to the flux tube. 



As shown in section I6.3[ the vertical extent of the flux tube (a = 190 km) is smaller 
than its lateral size (I = 360 km). Therefore, the cross section of the flux tube has a 
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flattened shape. iMagaral (1200 ll ) studied the emergence of large flux tubes with larger field 
strengths which form active regions. He found that the top of the loop decelerates when it 
reaches photospheric layers while the bottom of the flux tube continues to rise. The reason 
is that the photosphere is a convectively stable layer. Therefore the cross section of the tube 
changes from the circular shape to horizontally e xtended shape . The same process may 



Steiner et al. 



(120081 ) pointed out that an 



take place for granular-sized magnetic flux tubes, 
emerging flux tube can get pushed to the middle and upper photosphere by overshooting 
convection, forming atmospheric layers full of horizontal fields. This may also contribute to 
the flattening of the flux tube. 



7.3. Magnetic flux and energy 



The horizontal mag netic flux is estimated t o be 3.1 x 10 17 Mx near the apex of 



the flux tube. Recently 



Ishikawa &: Tsunetal (12009a ) reported that the occurrence rate 



of THM Fs is about 1. 1 x 10 10 km 2 s 1 in a quiet Sun region. Also, 



Tsuneta et al. 



(I2008bh and 



Ito et al. 



( 120091 ) have shown the presence of horizontal magnetic fields 



ubiquitously in the polar region, and that there is no differ ence on the intrinsic magnetic 



field strengths between the quiet Sun and the polar region. 



Litea (120091 ) suggested that the 



uniform fluct u ation s of longitudinal magnetogram signals over the whole disk observed by 



Harvey et al 



(120071 ) are the same phenomenon as the horizontal fields observed by Hinode. 
All these results together indicate that THMFs would have the same occurrence rate all 
over the solar surface. Combining these observations, we conclude that THMFs could carry 
a total magnetic flux of 1.8 x 10 25 Mx per day to the higher atmosphere what corresponds 
to a total of 7 x 10 28 Mx in one solar cycle (11 year). This estimation is based on the 
assumption that all THMFs reach higher atmospheric layers and that all of them carry 
the same amount of magnetic flux (given in this paper). Although our calculations are 
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base d only upon a single event, similar values have been found also by |jin et al.l (120091 ) 



and 



Martinez Gonzalez fe Bellot Rubiol ( 120091 ). The magnetic flux carried by THMFs is 



much larger than the total ma gnetic f 



corresponds to about 10 25 Mx (jHarvey 



ux in sunspot regions during one solar cycle which 



1993|). 



Ishikawa fe Tsunetal (j2009al ) have estimated the magnetic energy flux carried by 
THMFs, and found that the amount of energy flux is comparable to the total chromospheric 
and coronal energy loss, assuming that all THMFs reach above the photosphere. In order 
to perform a more precise estimation of the total magnetic flux and the total magnetic 
energy that THMFs could carry to higher atmospheric layers, we need to know what 
percentage of THMFs reach the chromosphere or higher layers in addition to the precise 
estimation of the magnetic flux of each THMF. Such statistical analysis was done by 



Martinez Gonzalez fe Bellot Rubiol ( 120091 ). who pointed out that 23% of the horizontal fields 



that emerge with clear bipolar footpoints in the photosphere reach the low chromosphere. 



7.4. Disappearance of THMFs 



We show that the THMF analyzed in this paper is rising through the line forming layer 
of the Fe I 630.2 nm lines and that it reaches higher atmospheric layers. The flux tube 
essentially maintains its integrity during the rising motion. At the beginning of this paper, 
we addressed three different possibilities of the THMFs disappearance. The observed event 
is clearly into the third of the suggested mechanisms, i.e., the THMF goes through the 
photosphere reaching the layers above the line forming region of Fe I 630.2 nm lines. So me 



THMFs do no t show any conspicuous footpoints in their lifetime (jlshikawa et al 



Jin et al. 



2008 



20091 ) . Such THMFs possibly harbor weaker magnetic fields than the equipartition 



field strength. Thus, the convective motion may destroy the flux tube before it reaches the 
upper photosphere, or it may force it to submerge below the photosphere. 



-24 - 



The authors want to express their sincere appreciation to Bellot Rubio for his 
pioneering work, and thank him for allowing us to use the SIRGAUS software. The authors 
are thankful for the encouragement by Bellot Rubio, Ruiz Cobo, Del Toro Iniesta, and 
Orozco Suarez. 

Hinode is a Japanese mission developed and launched by ISAS/JAXA, with NAOJ as 
a domestic partner and NASA and STFC (UK) as international partners. It is operated by 
these agencies in co-operation with ESA and NSC (Norway). The financial support from 
GA AS CR IAA 300030808 to Jan Jurcak is gratefully acknowledged. 



-25 - 



At=Os At=130s At=260s 




Fig. 1. — Temporal evolution of the continuum I c (top) and the integrated linear /circular 
polarization signal (middle /bottom). The integration range for the linear polarization signal 
is from —21.4 pm to 21.4 pm from the line center of the Stokes / spectra at each pixel. The 
integration for the circular polarization is carried out from —21.4 pm to —4.28 pm and from 
4.28 pm to 21.4 pm, taking the absolute values of the Stokes V spectra. The sign for the 
obtained circular polarization signal is determined to be the sign of the peak Stokes V profile 
at the blue wing. The time cadence is 130 s. The pixels to be analyzed in this paper are 
shown with red and green contours (see text). 
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Fig. 2. — Left : Correlation coefficient for the combination of the SOT/SP wavelength 
per pixel and the wavelength at the first pixel (offset). The contours show the correlation 
coefficient of 0.929, 0.931, 0.933, 0.935, and 0.937. The cross represents the position with 
the maximum correlation coefficient of 0.938, in which the wavelength per pixel and the 
wavelength at the first pixel are 2.14 pm pix -1 and 630.089 nm. Right : The SOT/SP 
mean Stokes / profi le (see text), which is wavelength-calibrated compared to the solar atlas 



([Beckers et al 



19761 ). Note that these profiles are normalized to the continuum intensity. 
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Fig. 3. — Observed (dashed) and the synthetic (solid) Stokes profiles for pixel^l from the 
first scan (top) and pixel#2 from the second scan (bottom). At the bottom of each panel, 
the difference between the observed and the synthetic profiles is shown in percent. 
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Fig. 4. — x 2 maps in the parameter space of the position and width of the flux tube for 
pixel^l (At = s, left) and pixel#2 (At = 130 s, right). The cross shows the position 
where the x 2 is minimum. The yellow line indicates Ax 2 = 13.7 (la) and the red ones 
Ax 2 = 21.3 (2a) and 29.8 (3a). Blue contours indicate the reduced x 2 values of 1.5 and 2.0, 
while black ones correspond to 3.0, 4.0, and 5.0. We define the reduced x 2 as X 2 /(3n — p), 
where n is the number of the wavelength points, and p the number of free parameters. 
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Fig. 5. — Various physical parameters as a function of the optical depth obtained for pixel#l 
(top) and pixel#2 (bottom). The gray lines show the parameters with the smallest x 2 value. 
The black areas are created by plotting the distributions of plasma parameters within the 
la confidence regions. The dashed lines enclose the line forming region of the Fe I 630.2 nm 
lines. 
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Fig. 6. — x 2 maps equivalent to those shown in Fig. H] for the first scan (At =0 s). Pixel 1-6 
corresponds to pixel^l. Lower left panel shows the integrated circular polarization signals 
for these pixels. The green and red lines represent the azimuth direction at logr = — 1. 




Fig. 7. — Same as Fig. El but for the second scan (At =130 s). Pixel 2-14 corresponds to 
pixel#2. 
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Fig. 8. — Top panels: Cross-sectional view of the magnetic field strength of the THMF 
for the first {left) and the second (right) scans. White lines show the inclination of the 
magnetic field vector. Middle panels: Cross sectional view of the Doppler velocity. The region 
with field strength above 150 G is contained with the solid lines. The shown distributions 
correspond to the pixels that are indicated by green lines in the azimuth maps from Figs. [6] 
and El For each pixel, the atmospheric parameters with the minimum \ 2 are plotted. The 
bottom two panels represent the normalized continuum intensity. 
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Fig. 9. — Representation of the vector field and the surface used for the definition of the 
horizontal magnetic flux. 
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